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Abstract. Recently there has been interest in extending ghost-free massive gravity, 
bi-gravity, and multi-gravity by including non-standard kinetic terms and matter cou¬ 
plings. We first review recent proposals for this class of extensions, emphasizing how 
modifications of the kinetic and potential structure of the graviton and modifications 
of the coupling to matter are related. We then generalize existing no-go arguments in 
the metric language to the vielbein language in second-order form. We give an ADM 
argument to show that the most promising extensions to the kinetic term and matter 
coupling contain a Boulware-Deser ghost. However, as recently emphasized, we may 
still be able to view these extensions as effective held theories below some cutoff scale. 
To address this possibility, we show that there is a decoupling limit where a ghost 
appears for a wide class of matter couplings and kinetic terms. In particular, we show 
that there is a decoupling limit where the linear effective vielbein matter coupling con¬ 
tains a ghost. Using the insight we gain from this decoupling limit analysis, we place 
an upper bound on the cutoff for the linear effective vielbein coupling. This result 
can be generalized to new kinetic interactions in the vielbein language in second-order 
form. Combined with recent results, this provides a strong uniqueness argument on 
the form of ghost-free massive gravity, bi-gravity, and multi-gravity. 
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1 Introduction 

In the past few years there has been a revival of interest in theories of Lorentz-invariant 
interacting massive spin-2 fields, for reviews see [1, 2], This class of theories includes 
ghost-free massive gravity [3, 4], bi-gravity [5], and multi-gravity [6, 7]. The theories 
generalize the Fierz-Pauli mass term [8] to the non-linear level in a way that avoids the 
Boulware-Deser (BD) ghost [9], for proofs see [10-17], see also [18] for a perspective 
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coming from BRST. These theories are of interest for cosmology [19-51] (as well as the 
mass varying [52-57], F(R) [58, 59], quasi-clilaton [60-62], and extended quasi-dilaton 
[63-68] extensions) and have been applied to studying condensed matter systems using 
AdS/CFT [69-73]. It is also an interesting fundamental question to ask how spin-2 
fields can interact consistently. 

It is often useful to consider these theories in the decoupling limit, using the 
Stiickelberg formalism introduced in massive gravity by [74, 75]. The decoupling limit 
serves multiple functions. First, it simplifies the theory by focusing on the inter¬ 
actions that arise at the lowest scale, which for ghost-free massive gravity occur at 
A 3 ~ (m 2 Mpi) 1//3 . Second, the decoupling limit can diagnose the presence or absence 
of a BD ghost. If the BD ghost is present it typically manifests itself as higher deriva¬ 
tives on the Stiickelberg fields [76, 77], giving rise to an Ostragradski instability [78]. 
As shown in [3], for ghost-free massive gravity after diagonalizing 1 , the scalar modes 
are described by Galileon interactions [79], which have second order equation of mo¬ 
tion, reflecting the absence of the BD ghost in the full theory. Finally, the decoupling 
limit simplifies the discussion of the Vainshtein mechanism [80], which is needed in 
order to avoid the van Dam-Veltman-Zakharov (vDVZ) discontinuity [81, 82], For a 
review of the Vainshtein mechanism see [83]. The Vainshtein mechanism is a strong 
coupling effect that occurs when the curvature R ~ m 2 , or in the decoupling limit 
when the Galileon operators at A 3 become large. For a proposal of a description of 
the strong coupling dynamics in Galileon theories at energies above A 3 , see [84], The 
complete decouping limit of ghost-free massive gravity including the vectors was de¬ 
rived in [85, 86 ], and the decoupling limit for multi-gravity was discussed in [87]. 

In this work we shall focus on bi-gravity (and the limiting case of massive gravity), 
which can be described in terms of two vielbeins e“ and /“. A natural question is 
whether ghost-free bi-gravity is complete, in the sense that it is the most general 
theory describing one massive and one massless spin-2 held that avoids the BD ghost 
and exhibits a Vainshtein mechanism. It is known that the potential interactions 
cannot be generalized further without introducing a BD ghost [3]. However, given that 
the potential term breaks diffeomorphism invariance, it is not a priori obvious that 
bi-gravity needs to maintain the same the kinetic structure and matter couplings from 
General Relativity (GR). 

Indeed both of these possiblities have been considered in the literature. Kinetic 
interactions that are not of the Einstein-Hilbert form were studied in [88-91], and 
matter couplings where the matter fields are coupled simultaneously to both vielbeins 
were considered in [5, 92? -102]. Such extensions may admit interesting cosmological 
solutions [103-108], and may be relevant for partially massless gravity [109-122], dark 
matter [123], and for electrically charged spin-2 fields [124] (which could be useful for 
holographic applications [125]). Furthermore, if such extensions that do not excite a 
BD ghost exist, then we should consider them as part of ghost-free bi-gravity. 

1 There is one interaction in massive gravity that cannot be diagonalized with a local field redefi¬ 

nition, however the equations of motion are still second order. 
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However, there are also several arguments that these extensions contain ghosts. 
In the metric language in four dimensions, the possibility of new ghost-free kinetic 
interactions in massive gravity was ruled out in [126]. Since the metric and vielbein 
formulations need not give the same theory when considering modifications of the 
kinetic structure or matter couplings [88, 91, 100], it is conceivable that there can be 
new kinetic interactions in the vielbein language. In [91], it was shown that in the 
vielbein language in first-order form that ghosts arise. In this work we extend these 
no-go arguments to the vielbein language in second-order form. 

Recently, there have also been arguments showing that modifications of the mat¬ 
ter coupling give rise to ghosts [94, 127, 128]. In fact, we will argue that these no-go 
statements are closely related. As we will discuss below, the uniqueness of the potential 
term itself places strong constraints on the matter coupling, and the no-go arguments 
for the kinetic terms also imply no-go statements for non-minimal matter couplings. 

Collectively, these arguments (including the arguments given in this paper against 
kinetic terms in the vielbein formulation in second-order form) show that there are no 
strictly ghost-free modifications of the kinetic term and matter coupling. However it 
may still be possible to view these couplings in an effective field theory sense below some 
cutoff. In particular, as discussed in [129], if there are ghost-free genuine interactions 
coming from the matter coupling or kinetic interaction, it may be possible for these 
interactions to become strongly coupled while still remaining within the regime of 
validity of the effective theory. This requires that these interactions come in at a scale 
A sx . that we can consistently take to be below the cutoff of the effective field theory 
(defined as the scale at which new physics enters), in other words so long as we can 
take A s . c . < A c o .. This logic was applied to the linear effective vielbein matter coupling 
proposed in [94] (see also related discussions for the matter coupling in [127, 128] and for 
new kinetic interactions in [89]), in which matter is minimally coupled to the vielbein 

v; = ae; + (3f;. (1.1) 

In [94], it was shown that this coupling is ghost-free perturbatively around the mini¬ 
superspace, and it was shown that there was a decoupling limit in which the theory 
was ghost-free. Thus treated as effective field theories, these models could potentially 
be of interest for cosmology. Following this logic, it is important not just to determine 
whether or not there is a ghost, but at what scale the ghost arises. In this work, 
we will consider a different decoupling limit which contains ghostly operators, and use 
the result of this analysis to place an upper bound on the cutoff of this matter coupling. 

Summary of Main Results. 

• By performing an Arnowitt-Deser-Misner (ADM) analysis [130] in the vielbein 
language, we show that the most promising kinetic and matter coupling exten¬ 
sions lead to a BD ghost. For the matter coupling, our argument gives a different 
perspective on how the constraint structure is modified than the analysis of [100]. 
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• For a generic matter coupling, there is a ghost already at the scale A 3 . This 
is consistent with the recent results of [127, 128], though we will use slightly 
different assumptions. 

• For the linear effective vielbein coupling, we find there is a decoupling limit, 
different from the one considered in [94], in which the following ghostly dimension 
nine operator appears 

0(TO)(8fc)T = (1.2) 

where is the matter stress energy tensor on a Minkowski background and 
= d^B v — d v B^ and where B M and 7 r are the vector and scalar Stiickelberg 
fields. This operator appears in both the metric and vielbein formulations. It is 
possible to take a different scaling limit because the interactions from the matter 
coupling vanish in the decoupling limit considered in [94], as we will show using 
a Galileon duality transformation [131, 132]. 

• Beyond the decoupling limit, the existence of this operator allows us to put an 
upper bound on the cutoff 

A c .o. < (ck/3 ) -1 / 5 (m 3 M | 1 ) 1 ^ 5 , for the linear effective vielbein. (1.3) 
It is likely this bound can be improved by considering higher dimension operators. 

• Using a similar decoupling limit analysis we identify ghostly operators for 11011 - 
Einstein-Hilbert kinetic interactions in the vielbein language in second-order 
form. In particular we consider the promising class of kinetic interactions pro¬ 
posed in [89] as well as the Jordan frame of the linear effective vielbein matter 
coupling. As we will discuss, both were shown to be ghostly in the decoupling 
limit analysis of [126]. 

Outline. In Section 2 we will review the proposals that have been made for non¬ 
standard kinetic interactions and matter couplings. In Section 3, we perform an ADM 
analysis to show that the most promising candidate modifications to the kinetic struc¬ 
ture and matter couplings contain a BD ghost. I 11 Section 4 we consider a general 
non-minimal matter coupling and will show using a standard decoupling limit analysis 
that the linear effective vielbein matter coupling is the unique choice that is ghost free 
at A 3 . In Section 5, we focus on the linear effective vielbein coupling, and use the new 
decoupling limit to place an upper bound on the cutoff of this theory. In Section 6 , we 
consider kinetic interactions in second-order form in the vielbein language, and show 
that these suffer from ghosts as well. 

Conventions. We will work with a mostly plus metric signature. We symmetrize and 
anti-synnnetrize tensors with unit weight, so for example = \{T^ V + T Vfl ). We 
use Greek letters / 1 , v, ■ ■ ■ to represent tangent space indices and a, 6, • • • to represent 
Local Lorentz indices. We denote the inverse of a vielbein e“ by e£. Whenever we 
work in the massive gravity limit, to emphasize that /“ is a fixed reference vielbein, 
we will write (/ re f)^- 
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2 Review of extensions to ghost-free bi-gravity and massive 
gravity 

2.1 Ghost-free bi-gravity 

In this work we will focus on bi-gravity, by which we mean a theory of two interacting 
spin-2 fields, with no other gravitational degrees of freedom. We will consider these 
as effective held theories with some cutoff scale A co .. We will discuss the cutoff in 
more detail below, but suffice to say that the cutoff is no larger than the Planck scale 
A c .o. Mpi. As a starting point, we will assume that an effective held theory for 
bi-gravity has the following properties 2 : 

1. The action is Lorentz invariant. 

2. There is a maximally symmetric vacuum, which we will take to Minkowski for 
this work. That is, g pu = af pu = br t ] fll/ is a vacuum solution for constants a and b. 

3. Perturbatively around the Alinkowski vacuum state, the theory describes a mass¬ 
less spin-2 held (with 2 degrees of freedom) and a massive spin-2 held (with 5 
degrees of freedom). 

4. We are able to trust the Minkowski vacuum solution within the regime of validity 
of the effective held theory. 

In order for the Minkowski vacuum to be within the regime of validity of the effec¬ 
tive held theory, all degrees of freedom should have non-vanishing kinetic terms about 
Alinkowski space. Otherwise those degrees of freedom are infinitely strongly coupled 
around Minkowski, and we cannot trust the Minkowski solution within the regime of 
validity of the effective theory. Requirement 4 is automatically satished if we require 
that non-perturbatively, using an ADM analysis, bi-gravity contains seven degrees of 
freedom. Then assuming we have satished Requirement 3, all the degrees of freedom 
have non-vanishing kinetic terms around Alinkowski. 

Ghost-free bi-gravity satishes all the requirements above. In the vielbein formal¬ 
ism, the action is given by [4-6] 

S = M 2 e S EH [e\ + M}S EH [e\ + m 2 M 2 ot S pot [e, f} + S matt [e, ^ m ], (2.1) 

where M 2 ot = (M e ~ 2 + Mj 2 )^ 1 , and where denotes a generic set of matter helds 
coupled covariantly to the vielbein e. Here we have defined the Einstein-Hilbert action 
for a vielbein e“ 

Seh [e] = \ J & A x |e| R[e\, (2.2) 

2 It may possible to weaken some of these assumptons, so this may be viewed as a starting point 
to fix ideas. 
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where R[e] is the Ricci curvature scalar built out of e. As written, this action is in 
second-order form, in the sense that the connection is considered a function of the 
vielbein, rather than being an independent held. Explicitly, the connection is given by 

< = el (0- b c -O c ‘ b - 0‘ c “) , (2.3) 

where O abc = e aa e b/3 d[ a e c ^y We also define the ghost-free potential interactions 

„ 4 n 4 

S P ot[e,f} = -- / J^c n £ aia2 a 3 a 4 ]Je ai JJ /% (2.4) 

n=0 i=l j=n -\-1 

where wedge products between forms are implied, e a f b = e a A f b . We choose the c n 
to ensure e“ = 5“ and /“ = 5“ is a vacuum solution, or equivalently that there are no 
tadpoles (which in these conventions amounts to ^2 n c n = ^nUCn = 0). We will also 
demand Y2 n n2c n = —4 so that around hat space the massive mode has a positive mass 
given by m. Otherwise we will leave the parameters arbitrary. 

Finally, we can add a matter sector which we denote as SAattfe, Vv] to denote a 
generic matter sector where the helds m are coupled to the vielbein e. We can also 
add a separate matter sector S m att [/, Xj\ where the helds Xj were minimally coupled 
only to /, but there are no direct interactions between m and Xj- 

In the rest of this section, we will work in the massive gravity limit, for which 

Mf —> oo, M e = Mpi (massive gravity limit). (2-5) 

This decouples the canonically normalized huctuations of /, so /“ is hxed to its back¬ 
ground value, /“ —$■ (/ref)))- If there is a matter sector coupled to /, this matter sector 
should be scaled as well. In this section we typically take the reference vielbein to be 

Hat (/ref)“ = 5“. 

2.2 Non-minimal matter couplings 

Since bi-gravity and massive gravity are theories with two metrics, it logically possible 
for a single matter sector to be simultaneously coupled to both vielbeins or metrics. 
The possibility was recognized by [5]. There it was pointed out that a generic coupling 
to two metrics would likely ruin the constraint structure preventing the existence of 
the BD ghost, and would also likely be ruled out by tests of the equivalence principle. 
Nevertheless, there could still be special choices for the coupling to matter that could 
potentially preserve the constraint structure and lead to interesting phenomenology. 

In [? ], the duality of the bi-gravity action (2.1) under the interchange of the 
two vielbeins e -H- / was used to motivate the study of matter couplings that would 
also preserve this duality. In [101] it was argued that the couplings could be given a 
geometric interpretation in terms of a Finsler metric. Matter couplings to two vielbeins 
can also be motivated by dimensional deconstruction [133]. However, in [93], an ADM 
analysis was performed which showed that generically when a matter sector is coupled 
to two sectors, a BD ghost is present. This result was also obtained using different 
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methods in [94, 97]. 


There is a simple argument showing that we expect a general matter coupling 
obeying a weak equivalence principle to contain a ghost (see also [134] for a version 
of this argument). 3 The most general matter coupling obeying the weak equivalence 
principle is given by 

S = M 2 e S EH [e} + M 2 f S EH [f} + m 2 M 2 ot Sp 0 t[e, f] + S matt [V{e, f),1> m ] (2.6) 

where V (e, /) is an arbitrary, non-linear combination of the dynamical vielbeins e and 
/. Inverting the relationship between e and V and sending e —>■ e(V, /), we can write 
this as 

S = M 2 S EH [e(V, /)] + M 2 S EH [f] + rn 2 M 2 ot S pot [e(V, /),/] + S matt [V, V>m] (2.7) 

In this form, we see that we have modified the form of the potential interactions. As 
a result we expect the BD ghost to re-appear. 

In [94, 97] an interesting coupling was found that avoids this argument. This 
coupling does introduce new degrees of freedom [94, 96, 100, 102], However it was also 
shown in [94] that there is a decoupling limit in which no ghosts appear, and the ghost 
does not appear to quadratic order in perturbations around the mini-superspace. Thus, 
it was argued that this coupling could potentially be of interest phenomenologically. 

This coupling, which we refer to as the linear effective viclbein coupling, is given 
by minimally coupling matter to the effective vielbein 4 

v; = v; = ae ; + pf;. (2.8) 

This coupling evades the above argument because the ghost-free potential terms S'potfe, /] 
are polynomials of the vielbeins, and so under this field redefinition a ghost-free poten¬ 
tial is mapped to a ghost-free potential. If we impose the symmetric vielbein condition 
Vabe^fv = VabGuf^ then this coupling can be described in terms of an effective metric 

fjfu = a 2 g pu + 2 apg^y/g-'f „ + p> 2 f pv . (2.9) 

However, the symmetric vielbein condition can be modified by the presence of this 
coupling, as discussed in [99, 100], so the metric and vielbein formulations are not 
equivalent unless the symmetric vielbein condition is added as a constraint. 

The linear effective vielbein coupling can also be derived by demanding that, at 
one loop, matter loops generate a ghost-free potential interaction. As a connection 
between modifications of the matter coupling and kinetic terms, we note that, in fact, 

3 As we will dicuss in Section 4, giving up the weak equivalence principle does not help, because 
we lose conservation of the stress energy tensor. 

4 We use a lower case v to distinguish a linear combination from an arbitrary nonlinear combination, 
given by V 
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at one loop the linear effective vielbein matter coupling will also generate a new kinetic 
interaction 


IV l00 p D C \LMp } j d 4 x^^ + c 2 LM^J d 4 x^^R cS + 0((R cS ) 2 ), (2.10) 

where c\ and c 2 are dimensionless parameters and L = log(Mpi//r), where /1 is a renor¬ 
malization scale. The second term, \J—g eS R cS , is a new kinetic interaction that was 
not present at tree level, but that appears at the same order in the curvature expansion 
as the normal Einstein-Hilbert term. We also note, as discussed in [94], that the lin¬ 
ear effective vielbein matter coupling spoils the possibility of the potential term being 
technically natural, as can be seen from the fact that the potential term is renormalized 

~ Mpj. 

The requirement that, one loop, the matter interaction generates a ghost-free 
potential interaction was used by [98] to propose a generalized version of the linear 
efective vielbein coupling 


S) 1/4f/(e ’ /)J ’ (2 ’ n) 

where t/“ is a general (non-linear) combination of vielbeins. Based on the argument 
earlier in this section, we would expect this coupling to have a ghost because after a 
field redefinition we essentially have changed the form of the ghost-free potential. We 
will verify this below. 

As this work was being completed, in [127, 128] it was shown that the only matter 
coupling that does not give rise to ghosts in the decoupling limit is given by the effec¬ 
tive vielbein coupling. The results of this work are in complete agreement with what 
is found there, but is complementary in the sense that we make different assumptions 
about the matter couplings we consider (we do not assume a weak equivalence princi¬ 
ple). Furthermore, we will extend that work by considering a more different decoupling 
limit in which the interactions of the linear effective vielbein coupling do not vanish, 
and we put an upper bound on the cutoff scale for this coupling. 

2.3 New kinetic interactions 

In a related but separate development, the possibility of a ghost-free kinetic interaction 
for a massive graviton beside Einstein-Hilbert has been considered. We will consider a 
kinetic interaction to be a generalization of the Einstein-Hilbert term in that it contains 
two derivatives in unitary gauge, though it can be non-linear in e and /. Perturbatively 
around Minkowski space, the kinetc term for a massless spin-2 field is fixed (up 
to field redefinitions) to given by the Fierz-Pauli kinetic term [8] 

£ fp = 


( 2 . 12 ) 






In [88] (see also [135]), it was shown that, for a massive spin-2 field, in addition to 
the Fierz-Pauli kinetic term there are also ghost-free kinetic interactions rs_/ f) 2 II" for 
n > 2. In four dimensions there is one such interaction, of the form 

£ der = e^e»' v 'P'°'d il H vvl d ll ,H pp ,H (y<J ,. (2.13) 

The ghost-freedom can be verified, for example, in the ADM formalism. Once we couple 
H pv to matter, by standard arguments [136] we will inevitably be forced to add futher 
non-linear interactions for H. Thus it is important to ask if there is a non-linear version 
of Equation (2.13) that propagates the same number of degrees of freedom non-linearly. 

The first issue we face when modifying the kinetic term beyond the linear level is 
that we have to choose how to describe the spin-2 field non-linearly. We can identify 
several different possible approaches that can be taken in constructing an extension to 
the kinetic term: 

• Metric vs. Vielbein Language: Do we choose to write the theory in metric or 
vielbein variables? 

• First-order form vs. Second-order form: Do we treat the connection associated 
with the metric or vielbein as an independent field? 

For ghost-free massive gravity with an Einstein-Hilbert term and minimal matter cou¬ 
plings, the four different formulations are all equivalent. 5 For modified kinetic in¬ 
teractions, however, there are (at least) four different approaches to constructing a 
non-linear action, which may lead to different theories. 

In [126], a no-go theorem was given for massive gravity in the metric language 
in four spacetime dimensions in second-order form which ruled out the possiblity of 
ghost-free non-linear kinetic interactions besides the Einstein-Hilbert term. This argu¬ 
ment is reviewed in Appendix A. While the argument was given for massive gravity, the 
argument also applies to bi-gravity, as massive gravity is a limit Mf —)■ oo of bi-gravity. 
Thus this result shows that for bi-gravity theories, there are no kinetic interactions, at 
least in the metric language. 

Modifications of the kinetic term in first-order form in the vielbein language, 
partly motivated by dimensional deconstruction of the Gauss-Bonnet term [138] in five 
dimensions, were considered in [133]. There it was shown that new kinetic interactions 
give rise to new degrees of freedom in the decoupling limit. The Lorentz Stiickelberg 
fields, which are non-dynamical in massive gravity and bi-gravity [85], can become dy¬ 
namical once the Lorentz invariance of the Einstein-Hilbert structure is broken, leading 
to new degrees of freedom. 


5 Strictly speaking in the vielbein there are additional branches of the theory that are not equivalent 
to the metric language [137] , however in this work we will restrict ourselves to the branch of the theory 
where the BD ghost is not present. 
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In [89], it was argued that there could be kinetic interactions that are ghost-free 
in the decoupling limit, even if there is a ghost non-linearly. As an example, the 
linear effective vielbein matter coupling can also be described in Jordan frame, where 
it appears as a new kinetic interaction. Starting from 

S = M 2 e S EH [e ] + M 2 f S EH [f] + m 2 M 2 ot S pot [e, f] + S matt [v(e, /),</>m], (2.14) 

and inverting the relationship between e and v by sending e —>■ e(v, f) we arrive at the 
action 

S = [a _1 (u - j3f)\ + M 2 S EH [f} + m 2 M 2 ot S pot [v , /] + S matt [u, V’m], (2.15) 

where now the fundamental helds that we vary in the action are / and v rather than 
/ and e. The notation S pot [v,f ] indicates that the parameters of the potential term 
change under this field redefinition. It was argued that since the linear vielbein coupling 
is ghost free in the decoupling limit, this kinetic interaction should be as well. We will 
revisit this point in Section 6 , given the new decoupling limit we use for the matter 
coupling. 

Nevertheless, based on this argument, a generalized set of kinetic interactions was 
proposed in the vielbein formalism in second-order form that could also be ghost-free 
in the decoupling limit 

S = M 2 S eh [ e\ + M 2 S EH [ae+(3f] + M 2 S EH [f] + m 2 M 2 ot S pot [e, f] + S maU [e,‘^ m ]. (2.16) 

We note that this interaction can be interpreted as the one loop kinetic term generated 
by the linear effective vielbein coupling given in Equation (2.10). This interaction was 
considered for Dark Matter in [123]. 

In this work we will close the gap and show that in the vielbein language in 
second-order form, there is still no loophole to the metric argument. We will consider 
the kinetic term, and show that there is a decoupling limit where the ghost appears 
at A 3 . In fact, as we will discuss in Section 6 , this kinetic term was already shown to 
contain a ghost in this decoupling limit by the no-go analysis of [126]. 

2.4 The possibility of an effective field theory description 

As discussed in the beginning of this section, we take the point of view that bi-gravity 
is an effective field theory that breaks down at some cutoff scale A c o .. Of course in 
gravity we always imagine that the cutoff is no higher than Mp\, but in principle this 
scale could also be lower. For example, A co . could be the mass of new degrees of 
freedom that are necessary to make the theory consistent. From this perspective, it 
could be that asking that theory has five degrees of freedom at all energy scales or that 
the theory has five degrees of freedom doing a full ADM analysis gives too strong a 
requirement. 

In order for the theory to describe a Vainshtein mechanism, the cutoff should 
be larger than A 3 . This is because, as discussed in the introduction, the Vainshtein 
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mechanism relies on operators at the scale A 3 becoming large. So if we want our 
effective theory to have a Vainshtein mechanism within the regime of validity of the 
theory, we must have A co . > A 3 . In the language of [129], A 3 is a strong coupling scale, 
but is not the cutoff of the effective held theory. 

In ghost-free massive gravity and bi-gravity, ghost-freedom in the ADM sense 
means that it may be consistent to consider that the cutoff larger is than A 3 , since it 
is possible that full unitarity is maintained above this scale. 

However, for the new kinetic interactions and matter couplings, there typically 
are ghostly operators that arise at some scale A g h os t 6 - Since these operators break full 
(not just perturbative) unitarity, we can safely infer that there must be a cutoff for 
these theories A c o . < A ghost , at which the effective held theory description must break 
down. It is also possible for unitarity to be broken by some other mechanism other 
than the appearance of a ghost, however in this work will focus on ghosts. 

Therefore it is interesting to ask if it is possible for the new kinetic interactions 
or matter couplings to have a strong coupling regime, in which we could trust the 
strongly coupled interactions without being sensitive to an unknown UV completion. 
In the language of [129], we can ask if there are irrelevant but important interactions 
coming from the kinetic interaction or matter coupling. As a minimum requirement 
to be able to trust a given operator in a strong coupling regime, we should ask for 
that operator to not introduce a ghost. As we will see, even this minimal requirement 
is very restrictive. A generic modification of the matter coupling or kinetic structure 
breaks unitarity at the lowest possible scale because the pure scalar or scalar-matter 
interactions are ghostly. 

In [94], it was shown that the linear effective vielbein coupling is special, in that 
the pure scalar-matter interactions arising from the linear effective vielbein matter 
coupling gave rise to second order equations of motion for the helicity -0 mode of the 
graviton. As a result, it was argued that the scalar-matter interactions could be treated 
as irrelevant but important interactions coming from the linear effective vielbein mat¬ 
ter coupling. 

Thus, it was argued in [94, 127, 128] that we can potentially trust the predictions, 
including a strongly coupled regime, of an extended matter coupling (or, as argued in 
[89], of a new kinetic term) if we can work in a regime where A 3 A c o < A ghost . 
In this picture A co . is an unknown scale that is not predicted within the effective 
held theory, other than being smaller than A^ost'• From this perspective, it becomes 
important not just to determine whether or not one or more ghosts is present, but also 
to establish the lowest scale at which ghostly operators arise. 

6 More accurately they arrange at a range of scales A g h 0 st.n, we will focus here on the lowest scale. 

' This is similar to how the Higgs mass is not predicted with an effective field theory containing 
massive W and Z bosons, even though there is a calculable scale ~ TeV at which (perturbative) 
unitarity is broken in that effective theory. 
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3 ADM analysis in the constrained tri-gravity picture 


In this section, we consider the most promising proposals for extensions of the matter 
coupling and kinetic interactions in the vielbein language in second-order form. We 
show that these extensions introduce new d.o.f.s by performing an ADM analysis. 

While at first glance the new kinetic term in Equation (2.16) may look related to 
a bi-gravity theory, we will find it more useful to view it as a constrained tri-gravity 
theory. We can write the new kinetic interactions in four dimensions as 

S = M 2 S £H [e] + M]S E „[f] + M 2 S EH [t>] + m 2 M p 2 ot S pot [e, /] 

+ j d‘x A; (vl - cel ~ f>K) • (3.1) 

where are 16 Lagrange multiplier fields. Without the constraint, this theory would 
simply be bi-gravity plus a decoupled copy of GR. The presence of the constraint 
completely modifies the Hamiltonian structure, and reintroduces a ghost. 

Therefore, without the constraint, we have three copies of (diagonal) diffeomor- 
phisms (diffs) and Local Lorentz transformations (LLTs). The constraint breaks this 
down to a single copy of diffs and LLTs (if we allow A to transform appropriately). In 
the massive gravity limit Mf —> oo and / is fixed to its background value, then these 
diagonal copies are broken and there are no gauge symmetries. 

Taking the massive gravity limit (given in Equation (2.5)) of the action (3.1), we 
find a modified massive gravity action 


S = M 2 e S E „\e] + M 2 S E „[t>] + m 2 M 2 ot S pot [e, /„,] + 


d'xK (vl ~ < - Wr .(© ■ 


(3.2) 

We will directly perform an ADM analysis on the massive gravity limit and find that 
there is a ghost. This might be surprising because at first sight it might have appeared 
that the action in Equation (3.2) could arise as a limit of bi-gravity. However this is 
not possible. Instead we see that it is best viewed as a limit of a constrained tri-gravity 
theory. 

We can also view the matter coupling from this perspective 


5 = M 2 e S EH [e] + M 2 S EH [f\ + m 2 Ad 2 ot S pot [e, / ref ] + *Smatt [v > ^m] 
J d-zAC « - ael - P(U)l) . 


(3.3) 


We will outline the ADM analysis of the matter coupling in Section 3.2. The matter 
couplings were shown to have a ghost in the vielbein formulation in [100]. Our formal¬ 
ism gives a different perspective on how the constraint structure is modified, but the 
conclusion is the same. 


3.1 Hamiltonian analysis of new kinetic term 

In this section we focus on the massive gravity limit defined by Equation (2.5) where 
/ is fixed. It is easy to extend the analysis below to allow / to be dynamical. 
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A very useful consequence of writing things in the constrained form of Equation 
(3.1) is that the distinction between first and second-order form vanishes, provided 
that we do not subsitute the solution to the A constraints back into the action. To see 
this, note that the equations of motion for the connections related to e and v are just 
the normal torsion-free conditions. Working in first-order form simplifies the analysis 
in D = 3. 

3.1.1 D = 3 case 

Since the potential term is not crucial to the analysis, we will not include it (it does 
not change any of the counting arguments below). We will also set M e = M v = 1 in 
this section since it is not necessary to track these scales for the ADM argument. The 
action is 

S = f Sate WV + + J d 3 xA£ - ael - /J(/„,);). (3.4) 

where /“ is some fixed tensor that we do not vary, and the Riemann two form is defined 
by 

R[u] ab = du ab + io ac io c b . (3.5) 

Passing to the Hamiltonian we find 

n = e a 0 P a [co] + ojQ b M ab [e, u] + v a 0 P a [fi} + ^ b M ab [v , n\ + A£C>, e, / ref ], (3.6) 

where 

PM = £ abc £ ij R b j[uj} = £ abc £V (diuf + U bd uf) , 

M ab [eM £ ab c £ ij me,uj\ = £ ab (, d^ + ufef) , 

C a , = v; - ae; - ^(/ ref )“- (3.7) 

The momentum conjugate to e“ is 

e>: =► < = -£abc£ l3 Uj b - (3.8) 

Using the canonical Poisson bracket between e and the conjugate momentum tt we can 
derive the relationship 

(3.9) 

Thus we identify several non-dynamical equations 8 (that is, equations that come from 
varying with a field that does not have a kinetic term) 

P a [u\ -a\° a « 0, 

PailA + 0, 

Mab[e, U)\ ~ 0 , 

M ab [v,/j,} w 0, 

v; - ae“ - /3(/ ref )“ » 0. (3.10) 

8 We distinguish between ‘constraints,’ which are relationships between dynamical variables, and 
‘non-dynamical equations,’ which allow for the possibility that the equation can be solved for an 
auxiliary field like the lapse ejj instead of for a dynamical variable. 
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Taking the Poisson bracket of the above conditions with the Hamiltonian we obtain 
further non-dynamical equations 

2 ^ b P b [u] - Q£ «AMA“ « 0, 

2 ^,rt]+e y AMAJwO, 

-2 e l oP b] [u] - a\ [ *e? « 0, 

-2 uJ“P%]+a!\ 6] «0 , 

- a:-Dfc[w]eg - R b v\ + au o yfc « 0, (3.11) 

where Dfc[u;] is the Lorentz covariant derivative using the spin connection u. 

Not all of the C“ have generated further non-dynamical equations. The reason is 
that there is no momentum conjugate to efj, so Cq trivially commutes with the Hamil¬ 
tonian. 

It is now useful to perform a counting argument. In three dimensions, the phase 
space contains 45 fields 9 

9 A“ +6e“ + 6 ujf + 6 vf + 6 nf + 3 eg + 3 cR + 3 v a Q + 3 R = 45 fields. (3.12) 

Meanwhile, there are 39 non-dynamical equations. 

Then, assuming we can solve all of the non-dynamical equations for all of the 
auxiliary fields (we can do this perturbatively), we are left with a phase space 

dimension of phase space = 45 — 39 = 6. (3.13) 

A healthy massive graviton in D = 3 has 4 phase space degrees of freedom, so this 
counting is consistent with the existence of a BD ghost. 

In ghost-free bi-gravity, the equation for the lapse e{] is a constraint after inte¬ 
grating out the shift, and the auxiliary components of the vielbein [6]. This structure 
is broken by the A constraint, however. After integrating out all of the auxiliary fields, 
the equation for the lapse e[j can be solved for the lapse, rather than being a constraint. 
This occurs because the equation that determines fx involves the lapse and the shift in 
a highly non-linear way. We will demonstrate that this occurs in the mini-superspace 
below. 

If all of these equations can be solved for the auxiliary variables, then there 
are no constraints. Commuting these conditions with the Hamiltonian then generate 
equations for the time derivatives of the auxiliary variables, there are no constraints. 
In other words, if we can solve for the lapse 

N = N(e“,u?). (3.14) 

Then the Poisson bracket of N with the Hamiltonian will simply give us an equation 
for N, it will not give us an independent constraint on the dynamical variables e“ and 
ojf b . In other words, the Dirac procedure terminates with no constraint to remove the 
BD ghost. 

technically there are 54 fields, if we include the momenta conjugate to the A“. This makes sense 
because the phase space should be even dimensional. However, these momenta will be immediately 
removed with 9 secondary constraints, so we ignore these fields. 
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3.1.2 Analysis for D dimensions 


Review of Massive Gravity. 

First it is useful to review the counting for massive gravity in D dimensions in 
first-order form. We will consider the generic case where the BD ghost is present. 
There are 


D 2 e“ + 


D 2 {D 


X ) u Cib 

— ^ 


D 2 (D + 1) 


fields. 


(3.15) 


There are no first class constraints because diffs and LLTs are explicitly broken by the 
reference vielbein. 

There are D equations of motion coming from eg. When a BD ghost is present, 
these are solved for the lapse and shift and there is no constraint among these equations. 

Then there are equations coming from varying tug 6 . These equations are 

constraints. One interpretation is that these constraints guarantee that the theory is 
equivalent to a metric formulation (which ultimately requires that D ^ D ~ 1 "> components 
of the vielbein remain unphysical). If these equations fail to be constraints, there are 
additional d.o.f.s in the phase space that we expect to be ghostly (even if they are 
not ghostly, they correspond to d.o.f.s that are infinitely strongly coupled around the 
Minkowski vacuum). We further expect these primary constraints to generate 
secondary constraints. The time derivatives of the secondary constraints will result in 
D(yD ~ l ' > equations that can be solved for the themselves. 

Lastly, there are second class constraints that remove some redundant components 
of the spin connection. In D dimensions, there are D(D — 1) spatial components of 
the vielbein e“, all of which have kinetic terms. However there are D( ' D ^ 1 > spatial 
components of the spin connection t which means that some of the spin connection 
components are redundant and are not momenta conjugate to the vielbeius. There are 
\D(D — 1)(-D — 3) non-dynamical equations that can be solved for these redundant 
vielbein components. 

To summarize, for massive gravity, there are: 


D + ^D(D — 1) + i D(D 2 — D + 2) non — dynamical equations. 

(3.16) 

Thus the total phase space is 


phase space dimension = — - — —- fields — -D(D 2 — D + 2) n.d. equations 


r D 2 — D — 2 


= 2 


+ 1 


d.o.f.s. 


(3.17) 


This is the right amount for a massive graviton in D dimensions + 1 BD ghost. Ghost- 
free massive gravity contains two additional secondary constraints that remove the BD 
ghost. 


New Kinetic Interaction. 
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The new kinetic interaction counting is now relatively easy to perforin. Because 
of the broken cliffs, we expect the counting for both e and v separately to follow the 
generic massive gravity counting above before imposing the constraints associated with 
A. Namely we expect 

fields in e, v,uj, p — n.d. equations = 2 D(D — 1). (3.18) 

Next we impose the constraints from A“. There are D 2 primary constraints. However 
as we have seen above, only D(D — 1) of these lead to secondary constraints. Thus 

D 2 - D - 2 

-2- + 1 ' 

(3.19) 

In order to remove the BD ghost, we would need that, after integrating out all of 
the nondynamical fields (that is, A,/i, v and the unphysical parts of e, u>), there was a 
constraint equation for e and oj. However given that this does not occur for D = 3 as 
we have seen above, it cannot occur for D > 3. We will also show this explicitly in the 
mini-superspace. 


dimension of phase space = 2 D(D — 1 ) + D 2 — D 2 — D(D — 1) = 2 


3.1.3 Mini-superspace 

To illustrate these points, we demonstrate how the above procedure works in the mini¬ 
superspace. This argument will hold in D dimensions. 

In the mini-superspace we write X\\ — L and A* = A A*. We also define lapses 
ef] = N, v o = K and scale factors e* = aSj,v l j = and the spin connections cof = 
coSj, /if = fid'j (the other components vanish). We will also introduce a matter source 
given by a perfect fluid with vanishing pressure that is minimally coupled to a. Then 
the action becomes 


S 



a D ~\D - 2) (ua - co 2 N) + b D ~ 3 (D - 2) (jib - fi 2 K^j 


-\-L (A — aN — (3) + A {b — ae — j3) + £ ma tt 


(3.20) 


The Hamiltonian is 


T-L = P 2 N + P^K + L(K — aN — /3) + A (6 -aa-fi)- Na D ~ 3 p(a). (3.21) 

where the momenta P a and Pb are given by 

P a = a D ~ 2 (D - 2)u, P b = b D ~ 2 (D - 2)/i. (3.22) 

There are four non-dynamical equations that follow immediately 

Ci = K - aN - p « 0, 

C *2 = b — aa — f3 ~ 0, 

C 3 = P 2 -aL- a D ~ 3 p » 0, 

<7 4 = P 2 + L w 0 . (3.23) 


16 







(3.24) 


The time derivative of C 2 generates another non-dynamical equation 10 


C 5 = {C 2 , H} = KP b - aNP a » 0. 


Solving C 5 for P b and C\ for K, we find that 

_ aNP a _ aN 
b ~ K ~ aN + p a 


(3.25) 


The key point is that this is non-linear in the lapse. In GR (and in bi-gravity), the 
equation C 3 is a constraint on P a . However, it now becomes an equation for the lapse 


C* = P + aPt = 1 


a " 


N 


aN + j3 


P: - 


a D ~ 3 p 


(3.26) 


Unlike the ordinary Friedmann equation, this is an equation that determines N in 
terms of p, rather than an equation for the scale factor. 

To summarize, we have demonstrated how the procedure outlined in the previous 
sections can be carried out in the mini-superspace, showing the existence of a ghost. 


3.2 Constrained vielbein picture for linear effective vielbein coupling 

The linear effective vielbein matter coupling was shown to lead to ghosts in the vielbein 
language in [100]. We outline an alternative argument demonstrating out current 
method to see that the matter coupling leads to new degrees of freedom in D = 3 (and 
thus in higher dimensions as well. Our starting point is 

S = M 2 S EH [e] + m 2 Mp 0t S pot [e, / ref ] + S mSLt t[v, if> m \ + j d 3 x A£(w“ - ae“ - /3(/ re f)“). 

(3.27) 

As in the previous section, we will ignore the mass term (taking m —> 0 ) as it does 
not affect the conclusions that we will draw, and we will take M e = M v = 1. Then the 
action is 

s = S EH [e ] + S'matt [u, vg + j d 3 x - ae“ - /3(/ ref )“). (3.28) 

We obtain non-dynamical equations 

P a - a\° a « 0, 

M ab « 0, 

-^(/ ref )“ » 0 , 

T^vl + A^ « 0. (3.29) 

We can first solve the last equation for A^. Then the first non-dynamical equation 
P a ~ aAg will involve the matter sector in a highly non-trival way, and this equation 
can be solved for the lapse, so it is not a genuine constraint. 

10 Commuting C 3 and C4 with H just generates the condition A = 0. 
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The time derivatives of these generate further consistency conditions. Note that 
very conveniently [P rna tt, H] = 0 since v has no conjugate momentum, so that the 
matter does not contribute to the consistency conditions in this language. 

2 u% b P b - ae ij Di[tjj] A“ « 0, 

-2ej+'M - C++ 1 « 0. 

a (-D k [cj]e o + OiU^ b e b k ) « 0. (3.30) 

There are 

9 + 9 + 3e[j + 3cJo = 24 non — dynamical fields. (3.31) 

We expect these to be determined by the non-dynamical equations. In addition there 
are 

6 e“ + 6c of 1 = 12 dynamical fields. (3.32) 

This leads to a counting 

15 + 9 primary conditions + 12 secondary conditions = 27 + 9 = 36 conditions. (3.33) 

We expect that we can solve 24 of these for the 24 auxiliary variables. Then the 
remaining 3 are genuine constraints on the auxiliary variables. If we are optimistic, 
these 3 constraints will generate 3 secondary constraints, so there will be 6 constraints. 
That leads to 6 phase space d.o.f.s, corresponding to graviton plus a BD ghost in the 
best case scenario. 

4 Cutoff for general non-minimal matter couplings 

We have shown in the vielbein language how two of the most promising candidates for 
extensions of the kinetic term and matter coupling modify the constraint structure of 
ghost-free massive gravity so that the BD ghost returns. Based on the arguments from 
Section 2.4 however, it may still be interesting to consider these extensions as effective 
field theories below some cutoff. The next three sections will address this question for 
a range of matter couplings and kinetic terms, including the linear effective vielbein 
coupling proposed by [94], 

In this section we show that a generic modified matter coupling will introduce new 
degrees of freedom at the scale A 3 . We will consider a broad class of matter couplings, 
and we will find that the only coupling that is ghost-free in the decoupling limit for a 
generic matter sector is for the matter to be coupled to a single linear effective vielbein, 
as proposed by [94]. There are other possibilities if we allow for restrictions on the 
choice of matter sector and couple that sector in a specific way. For example we can 
have a scalar field whose kinetic term is coupled to a single linear effective vielbein and 
with an arbitrary number of potential terms coupled to different effective vielbeins 
(confirming the results of [93]). 
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We consdier a generic matter coupling which is built out of an arbitrary number 
of “effective vielbeins” V[e, /]“ (which we allow to be non-linear functions of e and /, 
but not of their derivatives) 

S = M 2 S EH [ e ] + M}S EH [f]+m 2 M^ t S pot [e,f] + J>L'2 t [V,<U (4.1) 

V 

We do not assume that tyi V ; ] represents a minimal coupling between V and ip m , 
and a given matter field may couple to multiple vielbeins V. In other words, we do 
not assume a weak equivalence principle. As an example of what we have in mind, 
consider a vector field A M coupled to two metrics 

s m „, = f d 4 x|e| (AgftgftF^FA + J d l x\f\ , (4.2) 

where F tll/ = <9^ A, — d u A fJ . We will consider this coupling explicitly below and show 
that in fact it contains a ghost in the decoupling limit. We note that there can be 
matter couplings not of form we are considering because within each term the matter 
may only couple to a single vielbein. For example our assumption does not include the 
quasidilaton [60] and extended quasidilaton extensions [63]. 

Our goal is not to perform an exhaustive search of all possible matter couplings. 
A search was performed in [127] assuming the weak equivalence principle, and ruled 
out any coupling that is ghost-free in the decoupling limit besides the linear effective 
vielbein coupling. Our goal is to identify the sources of ghosts in the decoupling limit. 
We will see how the known cases that are ghost-free evade these problems, and we will 
also see how several candidate extensions go wrong. However, we do not claim that we 
have found the most general matter coupling that is ghost free in the decoupling limit. 11 

We expect a generic matter sector to introduce new degrees of freedom even in 
the decoupling limit, based on the argument in Section 2.2. To reiterate, if the weak 
equivalence principle is satisfied we can always do a field redefinition to make the cou¬ 
pling to matter a standard one. For a general matter coupling, this field redefinition 
will dramatically alter the form of the potential term, and in that description it will be 
clear that a ghost is present in the decoupling limit. This problem only gets worse if 
the weak equivalence principle is violated. Nevertheless it is interesting to see exactly 
where the problem lies in the matter coupling. 

As we will show below, the vielbein and metric formalisms are equivalent in the 
decoupling limit for matter couplings. In the metric language we can quickly see why 
we would expect a problem to arise. In the decoupling limit that we consider in this 
section, we will see that the metric and vielbein languages are equivalent. We can 

11 There are stronger requirements if one further demands that the graviton potential remains ghost- 
free after taking into account one loop corrections. For example, a scalar field with a kinetic term 
coupled to one metric and potential term coupled to another metric is ghost-free classically, but at 
one loop will detune the graviton potential. 
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expand a generic effective metric as gjf„ = rj /w + ciII Mi , + C 2 Hj lu + ■ ■ ■ , where c \ t 2 are 
constants and where 14^ = AThen the matter Lagrangian will contain the 
term 

£matt T) dfjd Q 7T d°3 V 7T, (4.3) 

where is the matter stress energy tensor on Minkowski space. There are also similar 
terms ~ <9 4 7r 2 T that arise from evaluating T on the non-trivial background of i r. Only 
special choices of the coupling could allow for second order equations of motion, given 
the inevitable appearance of higher derivatives on n in the action. Furthermore, since 
in the decoupling limit only the hclicity-0 mode n is coupled to matter, not the tensors 
or vectors B there is no held redefinition that can remove these higher derivative 
terms. 

4.1 Equations of motion 

Our approach to taking is to first compute the equations of motion for the Stiickelberg 
holds exactly, then take the decoupling limit. We will use superscript ^ or subscript 
(y) to denote that a particular quantity is associated with the effective vielbein V“. 
We assume that V can be expanded around e = f 12 

OO 

' , ('-/); = EE c S([ ir i* < 44 ) 

n =0 k 

where the c[\l are dimensionless parameters and where ([H n ]fc)“ refers to all possible 
tensor contractions of H and e at order n in 14, where H is defined as 

H; = e“ - /;. (4.5) 

We define the covariant derivative D^ v \ which has the schematic form 

D {v) ~ d + T {v) + J V) . (4.6) 

where T^’ M is the Christoffel connection associated with V (which acts on tangent 
space indices) and c o^’ ab is the spin connection (which acts on Lorentz indices). We 
will see below that u^' ab only enters the calculation for convenience, so we will choose 
to introduce a connection such that = 0. 

We will introduce Stiickelberg fields following the conventions of [85] for later 
convenience 

e“ = 

ft /;(*) = (4.7) 

where the tilde denotes a quantity that depends on the Stiickelberg fields. We define 
the stress energy tensor with respect to V = V (e, /) as 

«V) = _ M IA 

^ ~ \V\ SgW’^ ’ ( ' 8) 

12 We assume that V is analytic in a neighborhood of the point e = /. 
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where \V\ = det(l^), and where 


«£? = ruKVi. 


The equation of motion for the Stiickelberg helds can be written 
SS[e, /, -01, 

- = rn ivi pox o^ pot f 

v 


$4> a 


= m 2 ^p 0t £J lP0 , + = °’ 


(4.9) 


(4.10) 


where £? pot is the usual contribution to the Stiickelberg equation of motion from the 
potential term, 


£*,pp, = \f\DjP 


where 


"/xz/ ra \ 

(4.11) 

i[e,/] 

■ 

(4.12) 


l/) I/I Sgl 

To compute we vary the matter action with respect to the Stuckelberg helds 

,«d2t[VV>„ 


= / dV 


s~g(n 


flU 


w v) 


flU 


d a x\V\Tjy'IV^D‘/W, 

OJh 


(4.13) 


where we have used the relationship 

Sf; = Dl f) sr = W + (4.14) 

After integrating by parts we can express the contribution of the matter sector associ¬ 
ated with V to the Stuckelberg equation of motion as 


qv> = -l/i I i*P 


~ sv c 

rp[l T fV M 

W ‘ d/ ?J 


+ 6 6 ) 1 {V),U V C 


(4.15) 


So far this is an exact expression. Note that u^T >,ab actually cancels out of this ex¬ 
pression, we include it only so that the diagonalized diff and local Lorentz invariance 
is manifest. 

4.2 Decoupling limit 

To take the decoupling limit we first perturb around Minkowski 

h a 

e a = S a -|- 

** ** 2 M e ’ 

k a 

f a = 5 a + M 

«/ U Lfc 1 


2 Mr' 


f 

B a e au d 7r 

+ -4U- + a7r 


mMp 0 t m 2 M pot ’ 
\ a 

_ e \ a b /mM pot _ 6 


mM x 


(4.16) 


pot 


- 21 - 















(4.17) 


We then take the decoupling limit by sending 

m —> 0, M e , Mf —> oo, keeping A 3 = (m 2 M pot ) 1//3 hxed. 

However, this does not fully £x the decoupling limit for the matter couplings. In order 
to take a decoupling limit, we need to choose a scaling for the parameters 2 that 
appear in the effective vielbein. In this section, following [85], we will choose to keep 
these parameters hxed in the decoupling limit. This is certainly a safe choice, as we 
are guaranteed that no interactions will be singular in this limit. Furthermore, this 
gives us insight into the interactions at A 3 of the full theory if these parameters are 
0 ( 1 ). 


4.2.1 Ghost for generic non-minimal matter couplings 

With this choice of scaling for 2, the decoupling limit in the matter sector amounts 
to the replacement 

va ->■ v (i, i + n )j) = ( in the matter sect ° r )> ( 4 - 18 ) 

n,k 


where ([n™]*,)^ refers to all possible contractions of n and 6 at order n, and where we 
have dehned 


n ^ = 


_ 


d u d u ii 


A| 


(4.19) 


The Lorentz Stiickelberg helds completely decouple from the matter sector in this limit. 
Thus the metric and vielbein formalisms are equivalent in the decoupling limit, for the 
matter coupling considered in Equation (4.1). 


There are several useful simplifying features that occur in the matter sector when 
we take this decoupling limit. First, since e“ = e" = we can identify local Lorentz 
and diff indices. Thus in this section we will raise and lower indices with rj A1I/ . This will 
require us to be more careful in distinguishing a vielbein from its inverse, thus we will 
use the notation to denote the inverse of V p , so 



V p 

v V 


e 


(4.20) 


Note that V p will be symmetric, since it is built out of symmetric matrices that com¬ 
mute. Another special feature is that \V\ is a scalar, not a tensor density. This is 
because non-linearly a ratio of densities |H|/|e| is a scalar, but in the decoupling limit 
limit |e| = 1. Finally, we note that many of the connections vanish. In particular, 
r[ / A /i — uj\f' ah = = 0. This last condition is very useful. With this the torsion- 

free condition for / becomes 



8j; + f w* ft + = aj; + = o. 


(4.21) 


Finally, we note that fit) " is a physical quantity in this limit (that is, it need not 
only appear in covariant derivatives or in curvatures), since a difference in connections 
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— is a tensor, and 


0. 


To compute the equation of motion for the helicity-0 mode 7r, we write the vari¬ 
ation of the full action as 


SS 


d 4 


x 


f a 

(-'(f), pot 


J2 s kv )) 


V 


leading to the equation of motion for tt 


£ = -d 

( ~'7T 



= £ 


7T,pot 


+ ^2 £ - 


■K,(V) 


= 0 . 


v 


Plugging in the explicit formula for we hnd 


(4.22) 


(4.23) 


£ 


ir,(V) 


= d n 


\V\d p 


± (V),v 


V 


-i 



-ySf 


+ \v\r 



v 


-i 


"SVl 

T' 


(4.24) 


In general Equation (4.24) leads to a very complicated equation of motion that will be 
higher order in derivatives on n in the decoupling limit. Note that df ~ d 3 7i. Here, 
there are several sources of derivatives of /: fourth derivatives on n can arise from 
second derivatives acting on 5V/5f, or from derivatives acting on the connection h 
Third derivatives on n can arise from derivatives of 6V/Sf or from T ( ' h Of course, 
these are precisely the higher derivative terms we would expect from Equation (4.3). 
Furthermore, as discused above, there is no held redefinition that will remove the in¬ 
teractions that give these equations of motion. 


As a specific example, let us consider one of the couplings proposed by [98] 


V a 

V 


f K \ e \ + I/I 

V I/I 



(4.25) 


This kind of coupling was analyzed in the decoupling limit in [128], however we do the 
analysis here because it is instructive to see explicitly how derivatives on SV/Sf can 
lead to ghosts. Note that by construction, 


\V\ = K\e\ + \f\. 


In the decoupling limit, V"’ can be written as 


where 


e; = h(< 5“ + n“), 


n = 


(k + | i + n | \ 1/4 

V |i + n| ) 


(4.26) 


(4.27) 


(4.28) 
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The terms that lead to fourth order equations of motion that are universal (in the 
sense that they do not depend on the specific choice of matter sector) are d 2 5V/5f and 
<9T. Explicitly, these are 


c(d 4 i r) 


|i + n|ffWa a 3, 


ft 1 
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4|i + n 


- ft - 7 / 4 
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-i 


V 


-1 


per 


ft 1 


K 


4|i + n 


E " 1 

n- 7 / 4 )d a d,y p A. (4.29) 


This will lead to fourth derivatives on n that cannot be cancelled for generic matter 
sectors. This explicitly confirms the argument given in the introduction, that generic 
non-minimal couplings to matter will create a ghost in the decoupling limit. 

As a check, note that when k —> 0, VT —> /“ and ft — > 1. When we set n = 0 
above, all the fourth derivative terms cancel. This is guaranteed because of the torsion 
free condition. 


We can guarantee that the partial derivatives on the functions SV/5f vanish by 
requiring that V is linearly related to e and / 

V‘ = a v e% + p v f‘. (4.30) 

Then the equation of motion in Equation (4.24) becomes 

£,,m = f>v\V\d„ (vi DjPTfovf) . (4.31) 

In fact there is still a potential problem, which we can see by focusing on the contri¬ 
bution 

3 MV\% P (b< v) T^) . (4.32) 

In the last line we have kept all the terms that can potentially have fourth of 7r, coming 
from dT ~ d^n. If Du ^T IJL/ contains a connection, then this term will generically be 
non-zero. 


As an example, let us consider the case where we couple a vector field kinetic 
term to e and a mass term for the vector field to / 

Sm„, = J dV|e| + J dV|/| • ( 4 -33) 

This corresponds to summing over V = {e, /}, where we choose to take a e = (3f = 1 
and otf = /3 e = 0. Then 

T$ = m) (a„A„ - ■ (4.34) 
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where A 2 = g'AA^A^. Note that contains connections 

= 2m f »(/) s<7) (aA',% - Ao 


(4.35) 


A tedious but straightforward calculation shows that we can write the terms of the 
equation of motion that are fourth order in derivatives on i r as 

cr/3 r _ -i T[i 

g[j ) A^A a dpdpd^dr'K. 


cO 4 -) _ _ m 2 

C n,(V) ~ m f 
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-1 
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-1 


(4.36) 

These come precisely from the derivatives of the connections that appear in D^T^. 

To avoid this problem, we have several options. 

1. We can minimally couple matter to the vielbein e (so that for all V, /3y = 0). 

2. We can require a weak equivalence principle, so that we only couple to a single V, 

in which case the equation of motion for the matter fields will be = 0. In 

other words we can minimally couple to /, or to a single linear effective vielbein 
coupling. This is consistent with the result of [127]. 

3. Alternatively, we can put restrictions on the matter field stress-energy tensor 
TpjP so that the higher derivatives in £ n ,(v) cancel. 

Let us discuss option 3 in more detail. The simplest way to eliminate the higher 
derivatives that generically arise in Equation (4.32) is to demand that T^P = 14 (l ] PPP 
where U <x ^ is a scalar function of the matter fields. This automatically eliminates the 
terms that are fourth order in derivatives since there are no connections that appear 
in DjpT^y,. There are still potentially third derivative terms that could enter from 

derivatives of V 


S^v) = Pv\V\dJdM {V) \V\ 
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D/3 v \V\dM (v) dJ\V\ 
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(4.37) 


However, these potentially dangerous terms vanish, because 
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(4.38) 


In fact, this can be generalized slightly. The matter equation of motion can generically 
be written 

£ Id Al’ty) = O' (4.39) 

v 

Thus, so long there is only one stress energy tensor, which for definiteness we call Tp ), 
that is not proportional to then we can use the matter equation of motion to 

eliminate the contribution of this Tp) 


| = _Y \y\f)(V) 


(4.40) 
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This eliminates the problematic terms coming from , for example a 

contribution that would arise from a kinetic term. 13 Thus, a scalar field with a kinetic 
term coupled to one effective vielbein, and potentials coupled to any number of effective 
vielbeins, is ghost-free in the decoupling limit 

s = J d 4 x\v\ (- 1 - d 4 *l V\U (v) (x), (4.41) 

with each V given by a linear effective vielbein as in Equation (4.30). Based on the 
arguments in the next section, however, the only coupling of this form that will be 
completely ghost-free in will be the case studied in [127], where a the kinetic term is 
coupled to e or to /, and there can be a potential coupled to e and a separate potential 
coupled to /. 

There may be other ways to eliminate the higher derivative terms in S n ^vy Nev¬ 
ertheless, we emphasize that these couplings are not generic, and finding a T 1 "' such 
that £-k,(v) is second order in derivatives is not sufficient to guarantee that this TJyi 
comes from an action principle. We note that based on Equation (4.18), in the decou- 
pling limit we are considering in this section, a different way to frame the question is 
how to couple matter to Galileons in a ghost-free way. 14 

To summarize, we have found that for a generic coupling to matter, there will 
be ghosts at A 3 . The only known exceptions are either to place special restrictions on 
the matter coupling, or to couple a generic matter sector to a single, linear effective 
vielbein. 

5 Cutoff for the linear effective vielbein coupling 

In the previous section, if we did not want to put special restrictions on the matter sec¬ 
tor, we were led by a decoupling limit analysis to focus on matter couplings minimally 
coupled to a single effective vielbein 

>Smatt ‘S'matt ["0m W] • (5-1) 

where 

< = < + /?/“. (5.2) 

In fact, as we will show, the interactions of the linear effective vielbein actually vanish 
in the decoupling limit we have been considering. As a result, in this section we 
will perform a more detailed analysis tailored to the special linear effective vielbein 
matter coupling. In what follows we consider a new decoupling limit in which the 

13 This result can also be derived using the Galileon duality, since we can always go to a duality 
frame where the kinetic term is coupled to e where it is manifestly not coupled to the Stiickelberg 
fields. 

14 For example, one could consider multi-field extensions of the Galileons [139-141] as inspiration 
for matter couplings. 
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ghost will appear. This new limit exists (without causing any interactions to go to 
infinity) because the interactions for the linear effective vielbein coupling vanish in the 
decoupling limit we have been considering. 

To reiterate, our perspective on the decoupling limit is that it is a scaling limit 
that we can take in which the theory simplifies. When we take the decoupling limit, 
we will choose to scale the parameters a and f3 of the matter coupling to force the 
matter coupling to contribute non-vanishing interactions at A 3 . When we do this, we 
will find a ghost arises in the decoupling limit. After we have found these dangerous 
operators using the decoupling limit analysis, we can then turn back to the matter 
coupling away from the decoupling limit (that is, for general a and (3), and place an 
upper bound on the cutoff of the theory using the ghostly operators we have identified. 

5.1 Vanishing interactions in the original decoupling limit 

An action for a generic matter sector minimally coupled to the linear effective 
(5.2) can be written 

S'matt = J d 4 x\/-g( v '>C[g ( - v) ,ilj\, 
where we covariantly couple -0 to g( v \ The effective metric is 

9$ = Vabvffl. (5.4) 

In the decoupling limit we have been considering, we can write 

= (<* + I3f (’W + ^n„„ + nl) ■ (5-5) 

It is convenient to take a + (3 = 1. This can be acheived as follows: First rescale the 
coordinates to get rid of the overall (a + (3) 2 in front of the metric. Then, inside the 
brackets, we redefine f3/ (a + (3) —$■ f3. Then the effective metric becomes 

«£? = + 2 /?n„„ + /? 2 n= rue (y + pn“) + /m£). (5.6) 

Now under the Galileon duality transformation [131, 132] 

& = x 1 * - -fU^TT, 

A! 

<Vr d^x), (5.7) 

where x and p(x) are the dual coordinates and dual field respectively, the metric 
becomes 

9$ -»■ Vi*- ( 5 - 8 ) 

This field redefinition removes all interactions between the Galileon and matter, in the 
limit we are considering. 


vielbein 

(5.3) 
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5.2 New decoupling limit analysis 

In spite of this fact, we can still construct a decoupling limit in which the matter in¬ 
teractions do not vanish by choosing a different scaling for a and (3. 

Before taking the decoupling limit, it is worth massaging the expression for the 
effective metric into a form where the decoupling limit is easy to take. We are not 
interested in the helicity-2 modes (because as we will see these interactions will arise 
at a higher scale), so we set 


= A a 

fl = d, (x a + <f>“) 8*. (5.9) 

Writing A a b = 8 b + A a b for convenience (we are not perturbing at this stage, we are 
just reparameterizing A), the effective metric is 

9$ = Vat, « + /?/“) (ae b v + /3ft) 

= (a + /d) 2 ^ ( x a + /?<&“) d u (x p + fi$ p ) rj a p 

+aP [X afl X^ a + d^ a X au ], (5.10) 

where 

V = 6“Tj ac X c d 8t , (5.11) 

and where we have used the relationship 

X^v + X Ufl + Xfj, a X u a = 0, (5-12) 

which follows from A a c A b d r] cd = r/ ab . Inspired by the results of the previous section, we 
perform a diff on the metric 

<A> jfj(i) = d ll Y°d„Y< l g!$(x{Y)), (5.13) 

where 

W = + (5.14) 

After doing this, the effective metric becomes 

= (« + PfVtLv + olP [X atl d u $ a + d^ a X av ], (5.15) 

where now all the fields and derivatives are evaluated in the new coordinate system. 
The last term with the Lorentz Stiickelberg fields is an obstruction to removing the 
matter interactions with a coordinate transformation. 15 Note that this term is propor¬ 
tional to ct/3, so that it is not present if matter is coupled covariantly to a single vielbein. 


15 If we had carried the analogous argument out in the metric language, we would find similar 
non-vanishing terms, ultimately due to the fact that d^B^, is not a symmetric matrix. 



We will find it convenient now to take a + [3 = 1 so that the matter is canonically 
normalized. To take the decoupling limit, we use the scalings 
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a 
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<J>“ 



g A“ b / (mMpot) 


<5? + 


\ a 
A b 

mM pot 


B a d a n 

mMp 0t + m 2 M pot ’ 


+ ■■■ 


(5.16) 


where A“ b is antisymmetric. Taking the limit m —> 0, M po t —> oo holding A 3 = 
(m 2 M P ot) 1//3 fixed, we find (working in units where A 3 = 1) 


9$ = Vn» + K l - P) [Kfid a d u 7: + d M <9“vrA Ql ,]. (5.17) 


The Lorentz Stiickelberg held A should be determined by its equation of motion. The 
exact equation for the Lorentz Stiickelberg fields is given by [91] 


dC dC 

Vacfftp.rb VbcJn C\ J!n ~ 0 ’ 


(5.18) 


‘S/; 

where C = Ceh + £-mass + ^matter- For the effective vielbein coupling, this reduces to 

[m 2 M^Q(e, f)g^ + af)\v\T^} (ej/i - = 0. (5.19) 


where Q(e, /) is a scalar function that depends on the parameters in the potential 
term, whose precise form does not matter for our analysis (see [85] for a more detailed 
discussion of the decoupling limit of the potential term). Taking the decoupling limit 
of this equation amounts to keeping the terms that are 0(l/m) on the left hand side. 
This yields 

2(1,1 + n) [2A^ - F, u - X Mal d a d u) \ - f3(l - P) (-T Ma |cPct)7r) = 0, (5.20) 

where = d bl B l/ — d u B M . We can solve this perturbatively for A 

= 2^ lu/ (5-21) 

where p is a generic name for B, 7r, and the matter fields. 

Up to 0(p 2 ), we find that g cS becomes 


C = ^ + W 2 + 0(f\ ( 5 . 22 ) 

In this form, we see that g( v > is not simply a diff of flat space. In fact, the Riemann 
tensor associated with this effective metric does not vanish, because of the contributions 
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of the last term. This correponds to the existence of the ghostly operator that arises 
in this decoupling limit (restoring the dimensions) 



(5.23) 


where is the Minkowski stress energy tensor of the matter. Thus the equation 
of motion for the vectors given in Equation (4.15) will have higher derivatives arising 
from the connection associated with this metric, that cannot be removed with a field 
redehniton. 

This scalar-vector ghost is different from the normal ghosts considered in massive 
gravity, which involve purely scalar modes. Related to this, we note that it is normally 
assumed that we can ignore the vectors for solar system tests, because they are not 
sourced by matter. With this matter coupling, the vectors will sourced by T[^ if there 
is a background for n. 

Notice that to 0(<p 2 ) in the decoupling limit, there is no difference between the 
metric and vielbein formalisms. In other words, to this order, the modification of 
the symmetric vielbein condition in Equation (5.20) does not enter into the effective 
metric. Indeed, one can derive Equation (5.22) by using the metric language. There 
will be differences between the vielbein and metric languages at higher order, but it 
is unnecessary for our purposes to calculate these corrections since the ghost already 
appears at 0(<p 2 ). Of course, considering non-minimal covariant couplings (say between 
matter and R[g ^]) will only make these problems worse. 

5.3 Upper bound for cutoff scale 

In the previous section we have considered a special, simplifying limit of the theory 
where the ghost arose immediately at the scale A 3 . This limit is a good approximation 
to the full dynamics of the theory when a, (3 a/^3 /m ~ (Mpot/m) 1 / 3 1. As 

a result, we can rule out the possibility of taking a, f3 to be this large based on the 
analysis of the previous section. 

Of course it is not necessary to take a : (3 to be this large. Nevertheless, we certainly 
expect that the problems that appear in this limit will still be present even for more 
general values of a and (3. To make this more concrete, let us consider coupling to a 
massless scalar field 16 



(5.24) 


Without taking a decoupling limit, let’s consider the different kinds of interactions 
that appear when we perturb this around a Afinkowski background. We will consider 


16 Though none of the conclusions we make about the cutoff scale depend on us choosing a massless 
scalar field. 
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the theory in the metric language for simplicity. Up to quartic order, the action is 
schematically given by (after scaling coordinates so that a + /3 = 1) 


s ~ Jd‘ X (d X y + (^ + ^) (sxr + ox) 2 

+ J ^ t dBdd * (dx)2 + y ^} dB)2(dx)2 + 0{v% (6 - 25) 

where h^ u and k lw are the fluctuations in g IJV and respectively. We can see that 
if a,/I ~ (Mpot/m) 1 ' 3 then the decoupling limit from the previous section is a good 
approximation and the ghost will apear at A 3 . However, even for general cq/3, from 
this quartic action we see that scattering processes such as yy —» Btt are dominated 
by the interaction that we found in the decoupling limit 




(5.26) 


Because of these interactions are ghostly, we conclude unitarity is broken at or below 
the scale 

~ L 7 yP- < 5 ^> 

As a result, we can place a bound A c . 0 . < A g h os t- Since this operator also appeared 
in the vielbein formalism, this bound applies to the matter coupling in the vielbein 
formalism as well. Taking m ~ H 0 (where H 0 ~ 10 -33 eV is the Hubble parameter), 
and taking M pot ~ M P1 , we find 


A c .o. < (a/?)“ 1/5 x 10% ~ (a/3)” 1/5 x 10“ 9 eV ~ (a/3)“ 1/5 (100 m) -1 . (5.28) 


It seems likely that this bound can be improved by considering operators at higher 
order, however exploring this question is beyond the scope of this work. To discuss 
physics at distances shorter than A“ 4 in this class of theories, a UV completion is 
needed. In other words, A c . 0 . is a genuine cutoff, it is just not a strong coupling scale. 


The operator we have found is closely related to the perturbative ADM analysis 
of [94], In that work, integrating out the shift to sixth order in perturbations the 
Hamiltonian was found to be non-linear in the lapse 

n 3 ®: W* < 6 - 29 > 

It was suggested that we could identify N ~ ir /(Mpim 2 ) above. If we make this 
identification, then we would expect an operator of the schematic form 

° ~ (Si * )2 d #2 ' (5 - 30) 
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We see that this operator arises at the same scale that we found in Equation (5.26). 
It is interesting that we have recovered the same scale by different means. 

In fact, the operator we have found using a fully four dimensional analysis reduces 
to a covariantized version of Eq. (5.30) in two spacetime dimensions, in full agreement 
with the perturbative ADM analysis done in [94] in two dimensions. To see this, note 
that in two dimensions we can write 


F, u = jV, (5.31) 

for some scalar /. Then the important parts of the action are 

S 3 / {-\ F ^ _ = / dV (-if - 

(5.32) 

Integrating out / we find 

SD / d2x (^) 2(n '“ £ °' n)2 ’ (5 ' 33) 

which is a covariantized version of Eq. (5.30). 

Even though miitarity is broken at the scale A g h ost , it is still possible to view this 
matter coupling perturbatively below the cutoff scale. Since our (conservative) upper 
bound is parametrically larger than A 3 there is still potentially a regime of energies 
E where A 3 < E < A c . 0 . where the Vainshtein mechanism could potentially operate 
but we would still be within the regime of the theory. In addition, the fact that a 
ghost does not appear in the mini-superspace means that we can potentially trust this 
background to be within the regime of validity of the theory. 

However, we expect a ghost to receive a kinetic term on backgrounds that have 
anisotropy (and so can excite the above operator at quadratic order in perturbations 
around that background), including spherically symmetric backgrounds needed for 
solar system tests. 

Furthermore, if we are willing to view this coupling in a purely perturbative 
effective held theory sense (with no strong coupling regime), we should include all 
operators consistent with the symmetries present, not just this matter coupling. For 
example, while the matter sector is covariant with respect to the effective metric g ^ 
at tree level, this symmetry is broken by the mass and kinetic terms for e and /. 
Therefore beyond one loop one would expect to generate operators that broke the 
covariant structure of the matter coupling. 

Said differently, we cannot trust the interactions from the matter coupling in¬ 
volving the vector modes in a strongly coupled regime. All of the (non-vanishing) 
vector-scalar-matter interactions will give rise to higher order equations of motion and 
are ghostly. Thus there is 110 possibility of a strongly coupled regime for these inter¬ 
actions, a UV completion must be provided at this scale that removes the ghost, and 
physics above this scale will therefore be sensitive to the UV completion. 
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To summarize the results of the previous two sections, we have shown that for a 
very wide class of matter interactions, a ghost arises in the decoupling limit where we 
scale a,/3 ~ 1, in complete agreement with the recent results of [127, 128]. For the 
effective vielbein coupling proposed in [94], we have found that in fact, the interactions 
between n and matter vanish in this limit, which we have shown by doing a Galileo n 
duality transformation. Nevertheless, there is still a ghost for these interactions that 
arises as soon as they interactions become non-trivial. This can be seen by considering 
a different decoupling limit where we scale a,/3 ~ 1 j\Jrn with a + f3 fixed. We can 
therefore place an upper bound on the cutoff for these theories around a Minkowski 
background. 

6 Cutoff for new kinetic terms 

With the results of the previous section, we can now quickly consider the possibility of 
new kinetic interactions in second-order form in the vielbein language. We will work in 
the massive gravity limit (Equation (2.5)) to simplify the presentation, but the same 
conclusions when working in bi-gravity. 

6.1 Jordan frame of the linear effective vielbein coupling 

As emphasized in the introduction, we can interpret the new matter coupling as a new 
kinetic interaction. The Jordan frame of the linear effective vielbein matter coupling 
is 

S = M 2 e S E H[v(e , /ref)] + m 2 M 2 ot S pot [e, / ref ] + SAatJe,-0m], (6.1) 

where the linear effective vielbein is 

< = + «(/„ f)“- (6.2) 

The parameters 7 , k are related to the parameters a, /3 given in the previous section 
by 7 = 1/a, k = /3/a. The parameter 7 changes the canonical normalization for h (the 
fluctuation of e), so it useful to redefine M pot = 7 M e . 

We will perturb the vielbeins in the sense 1 ' 

7 = «+v „) G+M 

(/,sf); = e* + *“)C (6.3) 

Note that we can choose \ a b so that h^ u = is symmetric. This yields an effective 

metric 

9$ = d, (+ K$ a ) d v (xf* + ^) t ^/3 

+7 hfw + -jhl„ + K 7 [e Q/ A$" + d^ a e av \, (6.4) 

17 Note that we need to keep track of the flucation in e when working with the kinetic term. We 
will work with the vielbein fluctuation /i“ which is related to the metric perturbation by h v ; e ib e in + 
lh 2 — h 

4 "'vielbein ~ '‘metric- 
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where 


(6.5) 


= V" + 2^'“' ^ A)- 

The kinetic term can be expanded in powers of h 

M 2 e S EH [v] = M 2 e S EH [v ]\ h= o + M 2 e j d 4 x V (G H ^\h= o) + 0(h 2 ), (6.6) 

where = R ( 0 — \g$R^ Both of these terms are non-zero because of the piece 
proportional to firy in the effective metric. Note that in the previous section, the 
pure scalar interactions vanish in the kinetic term. In this picture, we do not need a 
Galileon duality transformation to see this, the argument is simply that the effective 
metric given has vanishing curvature when we only include the scalars, as can be seen 
directly from Equation (5.6). These two contributions have the schematic form 

/ M 2 m 2 M 2 m 3 

d 4 x(n'y) 2 —(^- 3 — d 2 ( Xddn) 2 + (k'j) — ^ — hd 2 (Xdd7i). (6.7) 

A-3 ^3 

The Lorentz Stiickelberg equation of motion is 

[™X 2 ot S(e, S)g^ + K M>|S ( ” )j "'] (epj - e‘/„“) = 0. (6.8) 

As before, when we integrate out the Lorentz Stiickelberg fields, the metric and viel- 
bein formulations are equivalent to 0(p> 2 ). Thus, we find a cubic interaction in the 
decoupling limit arising from the second term above 

0(pdiv)(dA)h — J 6 F lia d° d v Ti{£ h) ,w , (6.9) 

where £ is the Lichnerowicz operator. This will lead to higher order equations of mo¬ 
tion. 


While we have derived this operator directly in the language of the kinetic term, 
it also can be derived simply by starting from the Einstein frame picture and doing 
held redefinition on h of the form h —> h — dBddi r. We can also compute a cutoff scale 
directly in Jordan frame. We consider the same scattering process as before, Bit — > \X- 
The contribution of the above operator to this process (which also includes the matter 
coupling, /M pot ) breaks unitarity at the same scale we found before. Of course 

this is not surprising, as S-matrix elements are invariant under held redehntions. 

In fact, this kinetic term was already shown to contain a ghost in the metric 
language in the decoupling limit analysis of [126]. While the details of the decoupling 
limit are slightly different, nevertheless one can check explicitly using the formalism of 
that paper that the pure scalar interactions vanish, but the vector-scalar interactions 
do not. Because the metric and the vielbein are equivalent to leading order in ip, this 
also rules out the interaction in the vielbein. 

The above analysis can be easily generalized to bi-gravity, however one has to be 
careful to hrst diagonalize the kinetic terms for the fluctuations of e and /, since there 
is kinetic mixing with the modified kinetic term. 
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6.2 Einstein-Hilbert of the linear effective vielbein 

We can also consider the new interaction proposed by [89], which as shown in Equation 
(2.10) arises at one loop from the effective vielbein coupling. The interaction is 

S = M 2 S EH [e ] + M 2 S EH [v] + MjS EH [f } + m 2 M 2 ot S pot [v] + S matt [e, - 0 m ], (6.10) 

The only effect of the new term, to the order we are working, is to change the normal¬ 
ization of the fluctuation for e. The pure scalar interactions cancel, again because when 
we limit ourselves to the scalars then has the form of a diff of Minkowski space, 
so its curvature vanishes. This fact makes this kinetic term a promising candidate for 
a new kinetic term, a generic modification would already contain a ghost in the pure 
scalar sector. 

Thus we need to include the vectors. However, S E n[e\ does not contribute to the 
Lorentz Stiickelberg equations of motion at all, and so the equations are the same as the 
previous section. Again to 0(</? 2 ), the effective vielbein is the same whether calculated 
in the vielbein or metric formulation. In fact, S E h[g] does not contribute any vector- 
scalar interactions at all (since it is diff invariant), and so this kinetic interaction has 
the same ghostly operator we found in the previous section. 

6.3 Simultaneously modifying the kinetic term and the matter coupling 

Finally for completness we consider simultaneously modifying the kinetic term and 
matter coupling with different effective vielbeins 

S = M 2 S EH [v(e, /)] + MfS EH [f] + m 2 M 2 ot S pot [e , /] + SAattKe, f),ip m ], (6.11) 

where v and u are two different linear effective vielbeins. By doing a held redefinition 
on e we can write this in a form where the fundamental variables are, say, u and / 

S = M 2 S EH [v(u, /)] + M 2 f S EH [f] + m 2 M 2 ot S pot [u , /] + S ms , tt [u, ip m \, (6.12) 

where v is linearly related to u and /, the hat denotes the fact that the parameters in 
the potential term change under this held redehniton (but it is still of the ghost-free 
form). Then by the arguments of this section, there is a ghost unless v(u, f ) = u. 

7 Discussion 

In this work we considered non-standard kinetic terms and matter couplings in mas¬ 
sive gravity and bi-gravity, working in the vielbein language in second-order form. The 
bottom line is that all known extensions of the matter coupling to a generic matter 
sector and kinetic term contain ghosts. Furthermore, the ghost or ghosts always ap¬ 
pear when we take a decoupling limit in which the extension contributes non-trivial 
interactions. Beyond the decoupling limit, this translates into an upper bound on the 
cutoff for these extensions. 
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We gave an ADM argument that the most promising possibilities for these exten¬ 
sions change the constraint structure so that the BD ghost returns. These arguments 
are complementary to the ones in [ 100 ] and give a different perspective on why the 
constraints are lost. 

Additionally we considered the scale at which ghostly operators arise in this class 
of extensions. We considered a very general matter coupling, not assuming a weak 
equivalence principle, and showed that the generic matter coupling will give rise to 
ghosts at the scale A 3 , arising from interactions between the matter fields and the 
scalar modes of the graviton. There are special choices of coupling that avoid this 
conclusion, the only known choice that does not put further restrictions on the matter 
sector is the linear effective vielbein coupling of [94]. 

Furthermore, we found using a Galileon duality transformation that the effective 
vielbein coupling contains no purely scalar-matter interactions. By considering a dif¬ 
ferent decoupling limit where the matter coupling gives rise to genuine interactions, 
we found that a ghost arises in the decoupling limit and involves vector-scalar-matter 
interactions. Because of the presence of the vector mode, the operator we find rep¬ 
resents a qualitatively different kind of ghost than the one that usually arises in an 
arbitrary potential term. Related to this, the matter coupling will source the vector 
modes if there is a background for the scalar, unlike the case in massive gravity. 

The decoupling limit analysis we performed is strictly only a good approxima¬ 
tion to the theory if the parameters in the matter coupling are of order (M pot /m) 1//3 . 
Nevertheless, the ghostly operators we identified in this analysis will still be present 
in the full theory even beyond this regime. This allowed us to put an upper bound on 
the scale at which the theory breaks full unitarity and requires a UV completion. We 
expect that this bound can be improved by considering higher order operators. 

This analysis strongly suggests that the linear effective vielbein matter coupling 
cannot have a strong coupling regime. This is because the equations of motion of the 
vector-scalar-matter interactions break unitarity: the equations of motion are higher 
order, and these interactions cannot be removed with a field redefinition since they 
give rise to a non-zero scattering amplitude. As a result, we can at best consider the 
matter coupling or kinetic interactions to quadratic order around certain special back¬ 
grounds, however working to higher order in perturbations we will inevitably be led 
to interactions that are beyond the regime of validity of the field theory. Thus, while 
the matter coupling could be used to describe perturbations around an exact FRW 
solution, it cannot be trusted to describe structure formation processes that involve 
going to higher order in perturbations. 

This analysis can potentially explain some of the results found previously. It 
seems likely that the fact that anisotropies are required to see the ghost is related to 
the fact that the ghost arises from vector-scalar interactions. However, the precise 
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nature of the relationship between these statements is subtle, because the decoupling 
limit around Minkowski space of [89] is the same, yet that theory has a ghost in the 
mini-superspace. It could be interesting to explore this connection further. 


These extensions can still be considered perturbatively below A c . 0 .. Furthermore, 
iu the case of the vielbeiu matter copuling, the ghost does not arise perturbatively 
around the mini-superspace. 

This result, combined with other recent results [91, 100, 126-128], places strong 
constraints on the allowed form of ghost-free massive gravity, bi-gravity, and multi¬ 
gravity. It is remarkable that we do not need to invoke diffeomorphism invariance to 
conclude that we need the Einstein-Hilbert kinetic term and a covariant coupling to 
matter to avoid the presence of ghost degrees of freedom. 

Acknowledgments 

AAM is supported by an NSF GRFP. I would like to thank Kurt Hinterbichler, Nick 
Ondo, Johannes Nollcr, Rachel Rosen, Andrew Tolley for useful discussions, Raquel 
Ribeiro for useful discussions and comments on the manuscript, and especially Claudia 
de Rham for many useful discussions, comments on the manuscript, and encourage¬ 
ment. 

Appendix 

A Review of metric language no-go 

In [126], a no-go theorem was given for new kinetic terms in four spacetime dimensions 
for massive gravity (with the reference metric fixed as . The interactions 

considered were of the form 


£der(gj,dg,d 2 g), 


(A.l) 


with two total derivatives. 

Let us take a slightly different approach to explaining the analysis performed in 
that paper. We can construct a covariant form for the derivative of g by noting 



(A.2) 


where A is the difference in connections 


= r[sk - r[/k 


(A.3) 


Then the most general action is 
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where again, each term in £ der must have a total of two derivatives. The scale A der is 
in principle arbitrary. We want to study the leading order interactions, so we will scale 
A der in such a way that the leading order interactions survive in the A 3 decoupling 
limit 18 . 

We then choose to introduce Stiickelberg fields. We can introduce them through 

g as 

5V -t fjp, = d^d^gap^x)). (A.5) 


Then dehning 


Hfju — g^ v ^hiui 


we can express the action as 

n 


(A.6) 
(A.7) 


with all indices contracted with the flat reference metric rg w . Note that, unlike the 
analysis of [3], when the Stiickelberg fields are introduced through g the helicity-2 
mode h fW only ever appears in H^ u . Then perturbing 0" = x a + (mM P1 ) -1 £> a + 
(m 2 Mpi) -1 d a n , we can write H )W as 19 


. h^ u d^d^TT 

^ Mpi mMpi m 2 Mpi 

/ d^B a d^d a Tt \ / ()J3 n d u d a 7T \ 
ymMpi m 2 Mp\) \mMp\ m 2 Mp\) 


(A.8) 


with indices raised and lowered with We have neglected the fact that h^ u depends 
on 0 (which corresponds to replacing h lw —> 9 At 0 a (9 i/ 0 /3 /i ai a(0) above). In the decoupling 
limit that we consider, this can be shown to be irrelevant, for more details see [126]. 


We then demand that the equation of motion of h MI/ , B^, and n be second order. 
This is referred to as Property 2 in [126]. Applied to the derivative interactions, this 
is an extremely restrictive condition. It demands that the hclicity-0 mode n not enter 
into the decoupling limit at all, since any term in the action with n comes with at least 
4 derivatives per held (d 2 hd 2 7i), so the equations of motion are automatically higher 
order if they are non-vanishing. 

• First note that the Einstein Hilbert term gives rise to a kinetic term 

£fp ~ hd 2 h. (A-9) 

To be ghost-free, this must take the form of the Fierz-Pauli kinetic term. 

18 This choice of scaling for A<j e r is just for convenience. The real point is to study the leading order 
interactions. 

19 Note that there were typos in early versions of [126], this is the correct version of the Stiickelberg 
prescription that was used in the calculation. 
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• The first interactions that can arise in the kinetc term are purely 7r interactions. 
We can focus on these interactions by choosing A der ~ A 3 , 

/ d 2 n\ nn 

'^der,(m 2 M P1 ) 1 /3 ~ Ag 3 2 ■ (A.10) 

These interactions automatically give higher order equations of motion if they 
are not total derivatives. Furthermore there is no hope of mixing equations of 
motion at this scale. The coefficients must be chosen to remove these interactions 
entirely. 

• Once we cancel these interactions, we choose A der ~ (mMp J 1 / 3 , yielding 

Aler~ A 3 h 3 2 J (3B) 2 d 2 ■ (A.ll) 

These interactions again must be cancelled, which fixes more coefficients. 

• The next interactions arise at the scale A der ~ (mMpj) 1 ^ 6 . These have the form 

^der^mM^i/e ~ ^|<9 2 h<9 2 (dB ) 2 , d 2 (3B) 3 (A. 12) 

This scaling was not considered in [126]. However, we recently performed an 
analysis, new for this work, demanding that the equations of motion were second 
order when scaling A der ~ (mMpJ 1 / 6 (working with the most general lagrangian 
of the form (A. 7) up to quartic order in H). The only kinetic interactions that 
survive this decoupling limit analysis are (1) the Einstein-Hilbert term, and (2) 
the interaction £ 2 ,4 ~ eed 2 H i proposed by [88] (which is a total derivative in 
four dimensions). In order for (2) to be represent a genuinely new kinetic in¬ 
teraction, one would still need to End a ghost-free non-linear completion, which 
we do not exist since the decoupling limit analysis showed there was no non¬ 
linear completion for the closely related £ 2 , 3 - However this point is moot in four 
dimensions since £ 2,4 is a total derivative. As a result, in four dimensions, the 
decoupling limit analysis is strong enough to show that the Einstein-Hilbert term 
is the unique ghost-free decoupling limit up to quartic order in H. We will now, 
however, return to reviewing the analysis of [126]. 

Using this method, the most general action of the form given in Equation (A. 4) 
that has second order equations of motion in the decoupling limit up to quartic order 
in H. This four parameter family of interaction includes the Einstein-Hilbert term, 
as well as a candidate non-linear completion of the term proposed in [88], extended 
to cubic order. By construction, all of these interactions are trivial in the decoupling 
limit taking A der ~ (mMpj) 1 / 3 , but differ beyond that. 

£ d er, g.f. in d .l. = fl £a + &l£&i + ^2£fo 2 + C 1 £cd (A.13) 
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where the Lagrangians are all quartic order in H. C a starts at quadratic order C a ~ 
d 2 ( H 2 + H 3 + H A ). Meanwhile C bl2 start at cubic order ~ d 2 ( H 3 + H 4 ), and C CA is 
purely quartic order rv./ d 2 H\ 

The remaining four parameter family was then restricted by a perturbative Hamil¬ 
tonian analysis. In the mini-superspace approximation, two of the terms appear 
quadratic in the lapse. The final parameter was eliminated by considering anisotropies, 
where it was shown that the shift picks up a kinetic term and becomes dynamical. After 
a perturbative ADM analysis: 

= a {C a — C bl — 4 Cb 2 ) = O-C eh ■ (A. 14) 

which is just the Einstein Hilbert term expanded to quartic order. 

It was then argued that there are no interactions at any order in perturbation 
theory. This follows from the seemingly innocuous statement that at leading order 
in the helicity decomposition, any non-linear completion must be ghost-free. In other 
words, the interaction must be ghost-free in the decoupling limit at leading order in 
perturbation theory, without any possibility of mixing orders. 


In [89], it was pointed out that held redefinitions are apparently not covered in 
the above analysis. For example, we can imagine rescalings of H 



HfUJ —y X H^y. 

(A.15) 

Then if we 

introduce Stiickelberg helds through H , we find 



XH,y = X^~ + 2AT (/w/) + X^^y “, 

IVlpi 

(A.16) 

where 

,p _ d,By df.dy'K 

^ - mM-pi m 2 Mpi' 

(A.17) 


Because this decomposition does not preserve non-linear diffs, there will be higher order 
equations of motion on the Stiickelberg fields. Of course, these higher order equations 
of motion do not lead to ghosts, since all we have done is perform a held redefinition 
on Einstein-Hilbert. To make this explicit, we can perform a held redefinition on h to 
remove all higher deritvatives from the equations of motion 

h/iv h t w + A(A — l)Mpi'k /iQ 'k jy 01 (A.18) 

then rescaling h jlv —* h^/ A, B M —$■ B^/X, and tt —* tt/X. 

However, there is subtlety if matter is coupled to g IJiV = Tj tw + H /IU in the frame 
where the kinetic term depends on the metric g' = + A H^ IU . The same held 

redehniton we have done on to remove the higher derivatives from the kinetic 
term in Equation (A. 18) will introduce Stiickelberg helds into the matter sector. To 
analyze the decoupling limit in [126], it was implicitly assumed the metric appearing in 
the kinetic interactions was the same metric appearing in the matter coupling. Then 
demanding that no Stiickelberg helds appear in the matter coupling hxed the freedom 
to do held redehntions on H^ u . Similarly, we can consider arbitary nonlinear held 
redefinitions of H /1U . This can be absorbed into a (highly nonlinear) redefinition of 
h^ v . As before, this will introduce Stiickelberg helds into the matter coupling. 
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